Determining the upper limit of Fee for the Y (4260) 
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By fitting the R values between 3.7 and 5.0 GeV measured by the BES collaboration, the upper limit of the 
electron width of the newly discovered resonance Y (4260) is determined to be 240 eV at 9096 C.L., together 
with the BABAR measurement on Pee - B(Y (4260) + atx J/w), this implies a large decay width of Y (4260) > 


mm J[i final states. 


1. Introduction 


Recently, in studying the initial state radi- 
ation events, ete” — w«rsum'm J[V (YISR: 
initial state radiation photon) with 233 fb-^! 
data collected around „ys = 10.58 GeV, 
BABAR collaboration observed an accumulation 
of events near 4.26 GeV/c? in the invariant- 
mass spectrum of rtm J/wv [1]. The fit to the 
mass distribution yields 125 + 23 events with 
a mass of 4259 + 8*2 MeV/c? and a width of 
88 + 23*5 MeV/c?. In addition, the following 
product is calculated 


T'(Y (4260) > e*e-) - B(Y (4260) 2 ntr Jf) 


= 5.5 + rss eV. 


Since the resonance is produced in initial state 
radiation from ete collision, its quantum num- 
ber JPO = 177. However, this new resonance 
seems rather different from the known charmo- 
nium states with JPC = 177 in the same mass 
scale, such as (4040), (4160), and 1»(4415). Be- 
ing well above the DD threshold, instead of de- 
caying predominantly into DC) D) final states, 
the Y (4260) shows strong coupling to mt a7 J/v 
final state. So this new resonance does not seem 
to be a usual charmonium state. The strange 
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properties exhibited by the Y (4260) have trig- 
gered many theoretical discussions [2]-[12]. 

One suggestion is that the Y (4260) is the first 
orbital excitation of a diquark-antidiquark state 
([cs][es]) [2,3]. By virtue of this scheme, the mass 
of such a state is estimated to be 4.28 GeV/c’, 
which is in good agreement with the observation. 
A crucial prediction of the scheme is that the 
Y (4260) decays predominantly into D, Ds. 

Another opinion favors hybrid explanation [4— 
7]. In the light of the lattice inspired flux-tube 
model, the calculation shows that the decays of 
hybrid meson are suppressed to pairs of ground 
state 1S conventional mesons [13,14]. This im- 
plies that decays of Y (4260) into DD, D,D,, 
and D*D* are suppressed whereas D*D and 
D:D, are small, and D** D, if above threshold, 
would dominate (P-wave charmonia are denoted 
by D**). So it is interesting to search for the 
possible decay of Y (4260) > D1(2420) D. 

The third interpretation we like to mention 
is provide by Ref. [8], which suggests that the 
Y (4260) is the second most massive state in char- 
monium family. The author ascribes the lack of 
Y (4260) in e*e^ — hadrons to the interference 
of S-D waves, and also estimates 


T(Y(4260) + ete~) ~ 0.2 — 0.35 keV . (2) 


Besides the above interpretations, there are 
other kinds of proposals. The lattice study in 
Ref. [9] suggests that the Y (4260) behaves like a 
D,D molecule. In Ref. [10], it is proposed that 


the new state might be a baryonium, containing 
charms, configured by the A.-A,. In Ref. [11], 
the Y (4260) is considered as p-x«1 molecule while 
in Ref. [12], the Y (4260) is considered as w-x« 
molecule. However, all aforementioned specula- 
tions need further experimental judgment. 

Most recently, CLEOc collected 13.2 pb^! data 
at ys = 4.26 GeV and investigated 16 decay 
modes with charmonium or light hadrons [15], 
and the channels with more than 3o statisti- 
cal significance are mtx J/V (11o), nOn? J/y 
(5.17), and K* K-J/ (3.70). No compelling 
evidence is found for any other decay modes 
for the Y (4260), neither for the (4040) and 
(4160) resonances [15]. These measurements 
disfavor the p-x.; molecular model [11], bary- 
onium model [10], and high charmonium state 
explanation [8]. So far as other survived spec- 
ulations are concerned, such as charmonium hy- 
brid [4-7], tetraqark model [2,3], Di D molecule 
suggestion, and w-Xe1 molecule explanation [12], 
further experimental studies are needed to make 
more definitive conclusions. 

Since the Y (4260) was observed in e*e^ an- 
nihilation, it is expected that it contributes to 
the total hadronic cross section in e* e^ annihila- 
tion (or R values in other words). The most re- 
cent such data on measurements are from the 
BES experiment [18,19], as shown in Fig. 1 for 
Ec.m (= Vs) from 3.7 to 5.0 GeV. If we look in 
detail within the range from 4.25 to 4.30 GeV (re- 
fer to the inset of Fig. 1), it seems there is a bump 
around 4.27 GeV. Has this structure connection 
with the Y (4260)? We try to answer this ques- 
tion in this Letter. As there are other resonances 
nearby, we shall fit the full spectrum between 3.7 
to 5.0 GeV in order to get the information on the 
Y (4260). 


2. Fit the v-resonances 


The R values [18,19] used in this analysis were 
measured with the Beijing Spectrometer (BESIT), 
which is a conventional solenoidal detector ex- 
pounded in Ref. [16]. In the analysis below, the 
R values are converted into cross section a (e* e") 
by multiplying the Born order cross section of 
ete — wtp. These cross sections are plot- 
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Figure 1. R values for Ee.m. from 3.7 to 5.0 GeV 
measured by the BES collaboration [18,19]. The 
inset shows the R values in the vicinity of the 
Y (4260). 


ted in Fig. 2. There are clear peaks of 7)(3770), 
(4040), (4160), and w(4415). The data points 
have been used to obtained the parameters of the 
w-resonances [17]. 

To acquire the resonance parameters, we could 
simply fit the data set with cross section formula, 
each resonance with a Breit-Wigner 


1201} Ti, 


(3) 
where Pee, Ch, and T; are the mass indepen- 
dent electronic, hadronic and total widths, re- 
spectively, for a vector resonance of mass M pro- 
duced in the head-on collision of et and e^. No- 
tice that lee < T; in our analysis, the approxi- 
mation I; % T+ is actually adopted hereinafter. 
The summation of all the cross sections within 
the range we are studying is 


ails) = 37640). (4 


where indices 1, 2, 3, and 4 denote four resonances 
i0 (3770), (4040), (4160), and (4415), respec- 
tively. 
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Figure 2. Total hadronic cross section in nb ob- 
tained as c(e* e^ — hadrons) = R - 86.85/s (s 
in GeV?) from R values in Refs. [18,19]. Three 
sets of fit results, Fit 1, Fit 2, and Fit 3 corre- 
spond to three cross section forms 01, c», and 
703 as described in text. Two interference curves 
have been moved downward by 5 nb for display 
purpose, the dashed line at —5 nb corresponds to 
Zero cross section in the fit. 


In the 4»-family resonance region, if we assume 
that all two-body D™) D(9 states are decay prod- 
ucts of resonance, and not produced directly in 
continuum, we could therefore treat resonance 
and continuum incoherently. Nevertheless, for 
the four wide resonances, they are close and are 
expected to have some same decay final states, 
there must be interference between any two of 
the resonances. Therefore the amplitudes corre- 
sponding to each resonance, with the following 


form 
44/1204 rie 
Ts) = ——— ee 2 (5) 
[s — (M3)?] + i MIT? 
have to be added coherently to give the total am- 
plitude that, once squared, will contain interfer- 


ences of the type XT7T;. If the resonances are 


quite broad, the interference effect can also dis- 
tort the resonance shape, the width might appear 
broader or narrower, and the position of the peak 
can be displaced as well. In this case, the total 
Cross section 


ex(s) = |$ T;(3)| , (6) 


where T;(s) is given in Eq. (5). 

So far as amplitude is concerned, in principle, 
there are presumably relative phases between dif- 
ferent amplitudes besides the phase due to com- 
plex Breit-Wigner formula itself. So a more com- 
prehensive total cross section would be the sum- 
mation of amplitudes together with an additional 
phase, viz. 


2 


aa(s) = [2 Tile] (7 


Since what we actually obtain is the module 
square of amplitudes, only three relative phases 
could be detected in practice. 

The standard chi-square estimator is con- 
structed as follows 


m > (a°? (sj) — anan (8) 
2 — (Noev(s)? 

where g°®P (sj) indicates the experimentally mea- 
sured cross section at the j-th energy point, while 
gt^* (s j) i$ the corresponding theoretical expecta- 
tion at this energy point, which is composed of 
two parts 


athe (sj) = g"**(s;) + a^9^(s;) , (9) 


where c^?" denotes the contribution from contin- 
uum. Since there is little evidence in the data for 
any substantial variation of the continuum back- 
ground within the studied energy region, we pa- 
rameterize the continuum cross section with a lin- 
ear function 


ey (s) = A+ B(V/s — 3.700), (s in GeV’), (10) 


as has been used in Ref. [17]. Here we consider 
the continuum contribution as the background 


for measurement of resonance parameters, that 
is O49(s) = a^?" (s). 

In Eq. (9) o7®* denotes the contribution from 
resonances. The fit results are displayed in Fig. 2 
where Fit 1, Fit 2, and Fit 3 correspond to the 
three cross section forms 01, 02, and 03, as ex- 
pressed in Eqs. (4), (6), and (7), respectively. 
Although the synthetic curves for three fits are 
almost the same, as we expected, the interfer- 
ence effect deforms each resonance significantly, 
according to Fig. 2, w(3770) and (4040) be- 
come narrower while 1/(4160) and (4415) be- 
come wider when inference effect included. There 
exist constructive interferences as well as destruc- 
tive ones, and the interference behaviors for am- 
plitudes with and without extra phases are also 
very distinct. The fit results indicate that Tee is 
very sensitive to the fit strategy, and the largest 
difference between various amalgamation strate- 
gies could reach 5096. 

Since 4/(3770), (4040), (4160), and (4415) 
have similar decay feature, that is, all decay dom- 
inantly to DC?) DO) final states, we prefer the 
synthetic scheme of amplitude with phase which 
takes into account all possible interactions be- 
tween resonances. So far as the Y (4260) is con- 
cerned, the study of BABAR collaboration [1] im- 
plies that this new state may not be a common 
charmonium state, whose decay feature is rather 
distinctive from other family members, and this 
also obtains support from recent CLEOc mea- 
surements [15]. Therefore we add the Y (4260) 
contribution incoherently to the total cross sec- 
tions in Eqs. (4), (6), or (7). In this scheme, the 
upper limit of the production of the Y (4260) in 
ete” annihilation is determined, as expounded in 
the following section. 


3. Determination of lee of the Y (4260) 


Various fits to the data tell us that with 
limited knowledge on the nature of the reso- 
nances and comparatively meager data, we could 
only determine the Pee of the Y (4260) by a 
scan method. Specifically, we fix the mass 
of the Y (4260) at 4.259 GeV/c? measured by 
BABAR collaboration [1], and scan over the T; € 
(20,180) MeV/c? (8 MeV/c? step) and Te. € 


(0,500) eV/c? (10 eV/c? step) parameter space. 
In order to avoid some grotesque fit results due 
to random effect of the Y (4260) on the reso- 
nances nearby, in the scan procedure, the lower 
and upper bounds of the masses (total widths) 
are fixed to be 4100 and 4220 MeV/c? (30 and 
250 MeV/c?) and 4350 and 4500 MeV/c? (30 and 
300 MeV/c?) for (4160) and v(4415), respec- 
tively. 

Fig. 3 shows the likelihood as a function of the 
T; and Pee of the Y (4260), with the likelihood 
function defined as 


EET) = fn : ETE TI 
1 Ij —T,)? 11 
ER E L ( ) 
v2n0r, Or, 
where 
. ; 1 f f 
E =e (pC). am 


Here the X? is determined from Eq. (8), and fw 
an arbitrary normalization factor. The Gaus- 
sian term in Eq. (11) indicates that the measured 
width of the Y (4260) [1] is considered assum- 
ing I distributed as a Gaussian with the mean 
value T; = 88 MeV/c? and the standard devia- 
tion op, = 23.8 MeV/c?. Summing £,(Ti,, T?) 
with respect to I7, we obtain the variation of Ly 
versus Pee as shown in Fig. 3(b). The integral of 
the likelihood curve gives the upper limit of Tee 
of the Y (4260) at 9096 confidence level (C.L.): 


TY (4260) < 240 eve? . (13) 

By virtue of Fig. 3(a), we sum up £,(Ti,, T7) 
with respect to T$, to 90% fraction of the total 
area then obtain the upper limit of prid at 
9096 C.L. for each 0j So we obtain the 
variation of the upper limit of Fee versus T; as 
shown in Fig. 4. All the points are almost in a 
straight line, which indicates the ratio of the two 
quantities, or the upper limit of the branching 
fraction of Y (4260) — ete~ does not depend on 
the total width. Take the slope of the straight 


line in Fig. 4, we get 
B(Y (4260) > ete) < 24 x 1078, 
at 90% C.L. 
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Figure 3. Relative likelihood (Lw) distribution 
with respect to (w.r.t.) Tee and T; (a); and to 
Pee with T; integrated (b). The “N.I.” in (b) 
indicates normalized integral value for Lw. 


4. Other possibilities 


Although it is reasonable to treat the Y (4260) 
as an independent resonance in the fit and all 
the V-resonances are synthesized by amplitudes 
with phases, one may argue that other possibili- 
ties could exist since there is limited information 
about the properties of the Y (4260) as well as the 
other w-resonances. Thereby it's better to take a 
variety of effects into consideration. 

First, there are totally three categories or six 
schemes for summation of all the resonances in 
this analysis: 


1. The cross section of the Y (4260) is added 
directly with the other four resonances 
which are summed through cross sections, 
amplitudes, or amplitudes together with 
relative phases, as we discussed in Sect. 2. 
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Figure 4. The variation of the upper limit of Tee 
at 90% C.L. with respect to T+. The slope of the 
solid line denotes the averaged ratio T../T; = 
2.4 x 10-6. 


2. The amplitude of the Y (4260) is synthe- 
sized with other four resonances which are 
combined through amplitudes, or ampli- 
tudes together with relative phases. 


3. The five amplitudes corresponding to reso- 
nances 15(3770), (4040), (4160), (4415), 
and (4260), with relative phases are 
summed up to give the final synthetic cross 
section. 


Second, we consider the effect of different back- 
ground shapes in the fit. In Sect. 2, we adopt the 
first order polynomial to depict the background, 
nevertheless, higher order polynomial can be used 
to delineate the background as well. 

In addition, we notice one background shape 
once adopted by DASP group, who tried to take 
into account the threshold effect of the charmed 
mesons [20], 


6 
Odasp — Oete-—yutpu- ` (4o + 5 A; 6? F°’), (14) 
j=l 
with 
1 


et s/(3.1 GeV )? ’ 9) 


where Ag describes the contribution from contin- 
uum; j ranging from 1 to 6 indicates the DD, 
DD*, D*D*, D,D,, D;D*, D* D* thresholds, re- 
spectively. A; are free parameters; f; is the ve- 
locity of the relevant particles; and F is a over- 
simplified form factor defined above. 

The fit results show that the effect of poly- 
nomial background is at the same level with or 
smaller than the DASP background, so as an es- 
timation, we adopt the linear and the DASP back- 
grounds as two typical cases for background de- 
scription. 

Last, we also consider the possible energy- 
dependent T; of the Y (4260). There are argu- 
ments in some references [3-6] that the Y (4260) 
may be due to threshold effect just as the 
(3770), so the T; of the Y (4260) could be en- 
ergy dependent as follows: 


DOW). =T; 60V - 2mp;s)x 
Diyos 1 + (rDp.«)* (16) 
1+ (rpp-+)? Pp: l 
Here T; =T;(M), W = Vs is center-of-mass en- 


ergy, r the classical interaction radium, and p the 
D?- momentum, reads explicitly as 


1 
Dp; = 3y W? — Amie i (17) 


and p is the D3* momentum at resonance peak, 
viz. 
1 


em 2 V MP omis. (08) 


Dp:* = Pp** 


So T; in amplitude expression (5) should be re- 
place by the form expressed in Eq. (16) for the 
Y (4260). 

In a word, considering different summation 
schemes, background shapes, energy-dependent 
forms of the total width, a series of results give 
the upper limit of Tee ranges from 20 to 470 eV/c? 
at 9096 C.L. 


5. Discussion 


Our study indicates that the upper limit of Tee 
of the Y (4260) is about 240 eV/c?, which is con- 
sistent with the estimation presented in Eq. (2). 


Utilizing our upper limit Tee = 240 eV/c? for 
the Y (4260), together with the relation of Eq. (1), 
we obtain the lower limit of branching fraction at 
90% C.L. to be 


B(Y (4260) 2 rtr Ji) > 1.896. 


This means that the partial width I (Y (4260) > 
"1*7 J/v) > 1.6 MeV/c? at 90% C.L., which 
is much larger than the corresponding par- 
tial widths of w' (89.1 keV/c?) [21] and y” 
(44.6 keV /c?) [21,22], respectively. 

CLEOc measured the cross section for 
a*+a—J/y channel to be 58 pb^! [15], which is 
consist with BABAR result, 50 pb! [1]. As 
a conservative estimation, we regard the central 
value calculated in Eq. (1) is the same for CLEOc 
and BABAR, but adopt the improved accuracy 
prvoided by CELO, then we can obtain the fol- 
lowing lower limits at 90% C.L. 


B(Y (4260) > atar-J/p) > 19%, 
B(Y (4260) 2 1°n°J/y) > 0.796, 
B(Y (4260) > K+K-J/b) > 0.2%, 


and 
B(Y (4260) > XJ/w) > 3.4% , 


where X denotes (the numbers of the observed 
events by CLEOc are presented in parathses) 
ata (37), 909 (8), K+ K- (3), (5), 9 (1), n (0), 
nta 29 (0), and (1). 

As we notice up till now only results on 
hidden-charm final states were reported about the 
Y (4260), measurements involving open-charm are 
anxiously awaited to confirm existing specula- 
tions or provide clues for further theoretical in- 
quiry. In addition, more accurate Pee is still 
needed for a better understanding the nature of 
the Y (4260). 
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